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Abstract. Excitation of solar- like oscillations is attributed to turbulent convection and takes 
place at the upper-most part of the outer convective zones. Amplitudes of these oscillations 
depend on the efficiency of the excitation processes as well as on the properties of turbulent 
convection. We present past and recent improvements on the modeling of those processes. We 
show how the mode amplitudes and mode line- widths can bring information about the turbulence 
in the specific cases of the Sun and a Cen A . 
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1. Introduction 



Solar-like oscillations have now been detected i n a dozen main-sequence sta rs as well as 
in some red giant stars (see the recent review bv lBedding fc Kieldsenl 20061 ). Excitation 



of such oscillations is ensured by turbulent convection at the upper-most part of the 
convective zones. From the measurement of the mode amplitude and line-width, it is 
possible to infer the power supplied to the mode by turbulent convection. As highlighted 
in the present proceedings, this in turn permits to probe the turbulent convection in 
stars. 

Such measurements have been available for the Sun for several years and give the possi- 



bility to test the various proposed theoretical models of mode excita tion (e.g. lGoldreich fc Keelev 



1977tlBalmforth1ll992l:ISamadi fc GouplfeoOlllchaplin et al.1l2005t ). We consider here the 



model of ISamadi fc Gqupill (|200ll . SG01 hereafter) with the improvements proposed by 



Belkacem et all (|2006U B06b hereafter) 

The theoretical model of SG01 requires a prescription for the time-correlation between 
the turbulent el ements. On the basis of a 3D simulation of the upper part of the solar 



convective zone, Samadi et al. ([2003a, SNS03 hereafter) have shown that the eddy time- 
correlation, at large scale, is better fitted by a Lorentzian function (LF hereafter) than a 
Gaussian function (GF hereafter). Furthermore, excitation rates of solar p modes com- 
puted on the basis of the model of SG01 with a LF better agree with the excitation rates 
inferred from the helioseismic observations by IChaplin et aL (1998) than when using a 



GF. The open question is whether or not this result remains valid for other stars. 

The theoretical model of stochastic excitation requires a prescription for the fourth 
order moments (FOM hereafter) involving the entropy fluctuations and the turbulent 
velocity. SG01 assume the quasi-normal approximation (QNA) which consists in split- 
ting the FOM into the product of two second-order moments (SOM hereafter). In the 
solar convective zone, this simple closure model is however significantly biased as ver- 
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ified by the recent studie s performed by iBelkacem et aL ( 2006a , B06a hereafter) and 
Kupka fc Robinson! ( 20061 . and this volume). Furthermore, as shown in B06b, the skew 



introduced by the QNA leads to an under-estimation of the solar p mode excitation rates. 

Closure models more sophisticated than the QNA can be used. Among those, the so- 
called two-scale mass flux model (TFM hereafter) improved by iGrvanik &: Hartmann 
(2002) takes the asymmetries in the medi um into account but i s only applicable for 
quasi-laminar flows. B06a have generalized I Grvanik &: Hartmann! ( 2002f )'s approach by 
taking the turbulent properties of the medium into account. In order t o apply this model 
in the solar case, they have introduced the plumes dynamics following iRieutord fc Zahn 
(HH). Furthermore, as shown in B06b, the calculations based on this new closure model 
increases the contribution of the Reynolds stress contribution to the excitation rates of 
the solar modes. When the additional contribution due to the entropy fluctuations is 
included, the new theoreti cal calculation s fit ra ther well the maximum in the solar mode 
excitation rates derived by lBaudin et al. ( 2005 ). All the results by B06b are summarized 
and discussed in the present proceedings (see also Belkacem et al., this volume). 

Apart from the Sun, mode amplitudes and mode line-widths have b een derived for 
few st ars. Among those stars, a Cen A provides the best available data. iBedding et al 



(2004) lave derived the mode amplitudes and the mode line-widths for a Cen A . These 
mode line-wid t hs hav e been recently updated by Kieldsen et al. ( 2005! ). More recently, 
Fletcher et all (l2006h have derived the mode line-widths from the WIRE observations 



of a Cen A (jschou fc Buzasil 200lh . From those new measurements we have updated 
the work performed in Isamadi et aL ( 2004 ) by deriving new constraints on the rate at 
which energy is supplied by unit time to the p modes of a Cen A . These constraints are 
compared here with new theoretical calculations performed in the manner of B06b. 

This report is organized as follows: in § [2] we briefly recall the model of SG01 with 
the modifications proposed by B06b. We also present the different ways the eddy-time 
correlation is described and the different closure models investigated by B06a. Applica- 
tions of this improved excitation model are shown in § [3] to the case of the Sun (§ 13. ip 
and of a Cen A (§ 13. 2[) . Finally § [4] and § [5] are dedicated to conclusion and discussion 
respectively. 



2. Modeling the stochastic excitation 

The model of stochastic excitation we consider basically is that of SG01 with the 
improvements proposed by B06b: two terms are expected to drive stochastically the 
p modes: one term corresponds to the Reynolds stress, the second one is an entropy source 
term that corresponds to the advection of entropy fluctuations by turbulent motions. The 
energy supply per unit time into each mode is given by SG01: 

P^jjK + Cl) (2.1) 

where C\ and C| are the turbulent Reynolds stress and entropy contributions respec- 
tively. Their expressions are for radial modes: 



C|= / dmpofr dTe~ W0T d s r(w\wl)[r,T] (2.2) 

JM J-oo J 

C%= dmal^ dre- iU0T / d 3 r ((ws t ) 1 (ws t ) 2 } [r,r] (2.3) 

JM PO J-oo J 

where po is the mean density, w the vertical component of the velocity, s t the entropy 
fluctuation due to turbulence, a s = d(p/ds) p , m the mass enclosed in a sphere with radius 
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r, I the mode inertia, f r {£,r,m) = y~~§7j > ^ r ^ ne ra dial component of the displacement 
eigenf unction, q r is a functi on that involves the first and second derivatives of £ r (see Eq. 
(9) of ISamadi et al.ll2003b| ), finally r and r are variables related to the time and space 
correlation lengths respectively. 

Quantities labeled with 1 and 2 are taken at the spatial and temporal positions [xo — 
§> — §] and [xo + §, S] respectively: they correspond to two-point correlation products. The 
fourth-order two-point correlation product involving the velocity, < w\w\ >, is related 
to the fourth-order one-point correlation product, < wf >, as summarized in § 12.21 We 
point out that the fourth-order one-point correlation product, < wf >, corresponds to a 
FOM in the usual sense of Reynolds averages. In turn, the one-point correlation product, 
i.e. the FOM, is modeled as described in § 12.11 



2.1. Closure models 



The QNA (see iLesieurl 119971 Chap VII-2) allows us to calculate rather easily the FOM 
of w in terms of a product of the SOM, that is: 



(w 4 ) = 3 {w 2 f 



(2.4) 



The QNA is strictly valid for normally distributed fluctuating quantities with zero mean. 
However, the upper-most part of the convection zone being a turbulent convective system 
composed of essentially two flows, the probability distribution function of the fluc tuations 
of th e vertical velocity and temperature do not follow a gaussian law (see e.g. iLesieur 
19971 B06a). Hence the use o f the QNA is not really valid. This was verified by B06a 
and IKupka fc Robinsonl (j2006l ) with the help of 3D simulations of the outer layer of the 
Sun: Figure 12.11 shows the ratio between the FOM < w 4 > derived from a 3D simulation 
as detailled in B06a and the FOM derived from different closure models. As shown in 
Figure [2~T1 the use of the QNA under estimates, in the quasi-adiabatic region, by « 50 % 
the FOM of the velocity derived from the 3D simulation. 



A more sophisticated closure model, the mass flux model (jAbdella fc McFarlandll997 . 
MFM hereafter) , takes the effects of updrafts and downdrafts on the correlation products 
into account. Indeed, the presence of two flows introduces a skew when averaging the fluc- 
tuating quantities, since averages of fluctuating quantites over each individual flows differ 
from averages over the total flow. The basic idea of the MFM is to split any averaged tur- 
bulent quan tity 4> into two parts, one asso ciated with the updrafts and the other with the 
downdrafts. iGrvanik fc Hartmanri ( 2002 ') have proposed their so-called two-scale mass- 
flux model (TFM hereafter) which extends the MFM by accounting for the fact that hot 
drafts (or cold ones, resp.) do not necessarily coincide with the updrafts (downdrafts, 
resp.) . In order to close the third order moments and the FOMs, IGrvanik fc Hartmann 
( 20021 ) adopt the simplifying approximation that < (jf 1 >=< <j> > n where < . > denotes 
ensemble spatial (in the horizontal plane) and time averages. H ence, this approxima- 
tion c onsists in neglecting the turbulence within the flow. Finally, IGrvanik fc Hartmann! 
(2002) derive analytical expressions for the thir d-order moments and FOMs involving 
the temperature and the velocity. For < w A 
the following expression: 



Grvanik fc Hartmann! (|2002f ) obtained 



< w 4 > = (1 + sl ) < w 2 > 2 



with the skewness, S w , given by: 

S w ~ 



< ur > 



1 - 2a 



< w > 3 / 2 ^ a (l - a) 



(2.5) 



(2.6) 
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Figure 1. Fourth-order moment 
(FOM) of the velocity, < w 4 >, 
as a function of depth z, normal- 
ized to the FOM derived from the 
3D simulation. In solid lines the 
FOM calculated according to the 
CMP model, i.e. with the help 
of Eq. 12.71 and with the skew- 
ness, S m , given by Eq. 12.81 the 
dashed line is computed accord- 
ing t o the iGrvanik fc Hartmannl 
(2002)'s TFM given by Eq.[2~7lwith 
S m given by with Eq. l2.6l and finally 
the dotted line is the QNA, Eq.[2~4l 



where a is the mean fractional area occupied by the updrafts in the horizontal plane. In 
the QNA limit, i.e. when the random quantities are distributed according to a normal 
distribution with zero mean, we have necessarily Sw = . Hence, in the QNA limit, 



Eq. l2.5l does n ot match E q.fXH Grva nik fe Hartmannl (|2Q02r ) proposed to modify Eg. 12.5 
as follows (see Kupka fc Robinson 20061 for a discussion of this step) : 



< W 4 > = 3(1 



l5 2 
3 bw 



<w 2 > 2 



(2.7) 



Figure 12.11 shows that the FOM based on Eq. 12.71 with S w given by Eq. 12.61 results in a 
negligible improvement with respect to the QNA. However, when S w is derived directly 
from the 3D simulation and plugged into Eq. 12. 7\ Eq. 12.71 is a very good evaluation of 
the FOM derived from a 3D s imuation of the outer layer of the Sun as verfied by B06a 
and lKupka fc Robinson! (|2006h . 

B06a have generalized the TFM by taking both the skew introduced by the presence 
of two flows and the effects of turbulence inside each flow into account. Accordingly, they 
have derived a more accurate expression for S w , that is: 



1 



< w 2 > 3 / 2 



((1 - o)(l - 5a) Sw 2 - 3 < w 2 > )Sw 



(2.8) 



where Sw is the difference between the average velocity of the up- and the downdrafts 
and < w 2 > is the SOM. Both a and < w 2 > are supposed to be known. Finally, in 
order to close the sy stem, Sw needs to be mo deled. This is undertaken in B06a by using 
the plume mod e l of iRieutord fc Zahnl (jl995l ). More details can be found in B06a or in 
Belk acem et al. I (l2006ch . 



2.2. Two-point and Eddy-time correlations 

As seen in Eq. 12.21 and Eq. 12.31 the model of stochastic excitation relies on a prescription 
for the two point correlation products involving the velocity fEq. l2.2l and Eq. l2.3p and the 
entropy (Eq. 12.3ft . B06b proposed to relate the two-point correlation products involving 
the velocity to the one-point correlation with the help of QNA as follows: 



2 2 
< WjW 2 > = 



< llf > 
< W 4 >QNA 



< w\wl >QNA 



(2.9) 



where < w A > / < w 4 >qna is the ratio betwen the FOM of w and the one given by the 
QNA (Eq. I2.4[) . < w\w\ >qna is the two-point correlation product given according to 
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the QNA (see lLesieuij|l997l Chap VII-2), that is: 

(w?iu!) QN A = 2 ( Wl w 2 ) 2 + (wf)(wl) . (2.10) 

Note that the second term in Eq. 12.101 does not contribute to the excitation. The FOM 
in Eq. l2.9l is here computed according to the closure models presented in S I2.1I When the 
CMP is adopted, from Eqs. 12. 9[ Eq. 12.71 and Eq. 12.41 one derives the following expression 
for < w\w\ > (Eq. 7 of B06b): 

< w\w\ >= (1 + ^S 2 W ) < wjw 2 2 >qna (2.11) 

where < w\w2 >qna is given by Eq. 12.101 

SG01 introduce 4>ij(k,uj) the time and space Fourier transform of < W1W2 >qna (the 
first term in Eq. I2.10|) . SG01 modeled <fiij(k,ui) as follows: 

fa&w) = ^£ - ^ with E(k,oj) = E(k)xk(w) (2-12) 

where E(k,Lu) is the turbulent kinetic energy spectrum, &y the Kroenecker symbol, E(k) 
the mean kinetic energy spectrum and the eddy-time correlation function. Sev- 

eral eddy-time correlation functions have been investigated by SNS03. Among those, we 
compare here the Gaussian function (GF hereafter) and the Lorentzian function (LF 
hereafter). Finally from Eqs. l2~9ll2~T2l B06b derive the CMP version of Eq. [2^] (that is 
Eq. 10 of B06b). Note that as in SG01 and B06b, the entropy source term (Eq. [23]) is 
still modelled on the basis of the QNA. 



3. Results 

3.1. The Sun 

The calculations of the excitation of the solar modes is performed as detailed in B06b: 
the radi al eigenfuncti o ns, £ r , eigenfrequencies, an d the mode inertia, /, are those 
used by ISamadi et alj (|2003bl ) and originate from an ID solar model built according to 
Gough's (1977) non-local formulation of the mixing-length theory. The spatial and time 
aver aged quantities are obtai ned from a 3D simulation of the solar surface computed 
with lStein fe Nordlundl (|l998l )'s code and with a grid of 125 x 125 x 82. The total kinetic 
energy, E(k), in Eq. 12.121 and its depth dependence are obtained directly from the 3D 
simulation. Finally, quantities related to the convection (the SOM of w and St, and a) 
and thermodynamic quantities (the mean pressure p$, po and a s ) are also obtained from 
the 3D simulation. The results of the calculations of the mode excitation rates (Eq |2.1[) 
are shown for the Sun in Figure 13. 2i top panel. 

3.2. a Cen A 

Calculations of mode excitation rates are performed for a Cen A in a similar way as 
for the Sun: the simulation has a grid of 125 x 125 x 82, an effective temperature of 
T e g = 5805 ± 20 K and a gravi ty of logq = 4.305 in agre ement with the constraints 



derived from interferometry fsee iMiglio fc Montalbanll2005f ). We compute a ID model 



with the CES AM code assuming standard physics and the same T c ff and log g as for the 
simulation. Finally, £ r , u>q and / are computed with an adiabatic code from the ID model. 
The results of the calculations of the mode excitation rates are s hown for a Cen A in 
Figure |3~!?1 bottom panel. The seismic constraints are derived as in Samadi et alj ( 2004 ) 



except that, after adding the apparent amplitudes for t to those for £ + 2 (with £ — 
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Figure 2. Top: Rates at which 
energy is injected into the so- 
lar modes. The filled circles cor- 
respond to the heli oseismic con- 
strain ts obtained by iBaudin et al.l 
(2005). The lines correspond to dif- 
ferent theoretical calculations: the 
solid line uses the Lorentzian func- 
tion (LF) and the CMP, the dashed 
line is as the solid line with only 
the contribution of the Reynolds 
stress, the dot-dashed line uses the 
LF and the QNA closure model, 
the triple-dot-dashed line uses the 
Gaussian function (GF) and the 
CMP. Bottom: Same as top for 
a Cen A . The thick lines corre- 
spond to the constraints obtained 
fro m the observed spectru m derived 
by iKieldsen etHl 1)20051) and the 
av eraged mode line-widt hs derived 
by IKieldsen et all (120051. the thick 
solid line) and by iFletcher et akl 
(2006, the thick dashed line). The 
thick dotted lines correspond to ex- 
treme values derived from the esti- 
mated errors on the amplitude and 
line- width measurements. The thin 
lines have the same meaning as in 
the top panel. 



and 1), we divide the summed amplitudes by the sum of the visibility factors of all the 
summed modes including the multip lets. Furthermore, the mode masses (Ai = ) are 



evaluated as in lSamadi et al.l (|2004f ) at the photosphere 



4. Conclusion 



Comparison between recent helioseismic constraints obtained bv lBaudin et alj (|2005D 
and the theoretical calculations based on the quasi-normal approximation (QNA) con- 
firms the result of SNS03 that the Lorentzian function (LF) results in a better agreement 
with the observations than the Gaussian function (GF). 

As for the Sun, we find that for a Cen A the theoretical calculations based on a 
Lorentzian eddy-time correlation function is closer to the observations than those based 
on a Gaussian one. 

In the case of the Sun, B06b have shown that the QNA results in a significant under- 
estimation of the contribution of Reynolds stress to the mode excitation. On the other 
hand, the so-called closure model with plumes (CMP) proposed by B06a increases sig- 
nificantly this contribution. The theoretical calculations limited to the Reynolds stress 
only underest i mates by ~ 15 % the maximum in the excitation rates P derived by 
IBaudin et al. (|2005l ). On the other hand, when the entropy contribution is added, the 



theoretical calculations based on this improved closure model fit the helioseismic data. 
In the case of a Cen A , the theoretical calculations underestimate the observations by 
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a factor ten. We note also that the accuracy at which mode amplitudes and line-widths 
are measured is more or less of the same order as the difference between theoretical 
calculations using the QNA and those using the CMP. This emphasises the need for 
more accurate seismic data. 

Whatever the assumptions about the eddy-time correlation or the closure model, the 
theoretic al calculations are fou nd to be closer to the cons traints derived using t he line- 
widths of Fletcher et al. I (|20Q6l) than to those obtained for iKieldsen et all (l2005h 's. 

The space based mission CoRoT (see the recent review bv lMichel et al.ll2006h is the 
only asteroseismology mission that in the very near future will enable us to derive, for 
a large set of solar-like oscillating stars with different effective temperature and gravity, 
the rates at which energy is supplied to the modes by turbulent convection. The quality 
of the data is expected to be significantly higher than current observations. Hence, these 
data will thus very likely permit to discriminate between the best description for the 
eddy-time correlation function and closure models. 



5. Discussion 

For the Sun, the largest discrepancies between theoretical calculations and the obser- 
vations are seen at high frequency (y > 4 mHz) and at low frequency (y < 3 mHz). 

Part of the remaining discrepancies with the helioseismic data can be attributed to 
the closure models considered here which are not able to reproduce correctly the FOM 
from the 3D simulation in the super-adiabatic region. Indeed, as seen in Figure 12.11 all 
the closure models considered here overestimate the FOM. B06b have shown that the 
main effect on the excitation rates is at rather high frequency and results in an increase 
of the order of ~ 15 % when the CMP is adopted. 

In B06b, the one-point correlation of the velocity was generalized to a two-point cor- 
relation product using the QNA (see § 12. 2 [I . It was shown in B06b that, at small scale 
length, this approximation reproduces satisfactorily the constraints derived from the 3D 
simulation. At large scale length, however, this approximation overestimates those con- 
straints. The effect on the excitation rates of this overestimation at large scale remains 
to be evaluated. 

Finally, we must point out that the calculation of the mode excitation rates requires 
the calculation of the mode eigenfunctions. The latter are obtained on the basis of global 
ID stellar models. In such stellar models, the structure of the outer layer is very poorly 
modeled. Therefore, this poor description must have some (currently unknown) conse- 
quences on the inferred properties of the mode eigenfunctions in particular on the high 
frequency modes which are almost confined near the photosphere. In turn, as most of the 
excitation occurs near the photosphere, this poor description is likely to introduce some 
(unknown) biases on the predicted excitation rates in particular at high frequency. 

The discrepancies at low frequency can probably be attributed to the adopted model 
for the eddy time-correlation which is probably too simple. Indeed, the departures from 
a Gaussian behavior were attributed in SNS03 to the plumes, which are more turbulent 
than the granules. Near the photosphere, the excitation of the p modes is almost solely 
due to the plumes. However, one Mm deeper, the excitation is no longer dominated by 
the plumes. As the granules are less turbulent than the plumes, their time-correlations 
follow more likely a GF distribution than a LF. Then one Mm deeper the eddy-time 
correlation is expected to lie between the LF and the GF. This however remains to be 
checked. 

For a Cen A , the large discrepancies between the observations and the theoretical 
calculations remain to be understood. In addition to what is mentioned above for the 
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case of the Sun, the ways the mode excitation rates were computed in the particular case 
of a Cen A may have some (unknown) effects on the computed excitation rates. Indeed, 
the adiabatic eigenfunctions were obtained from a ID model based on the local mixing- 
length theory. Finally, the 3D simulation considered here uses solar abundances. On 
the observational side, we shall mention that the excitation rates were derived from the 
observations using mode masses evaluated arbitrarily at the photosphere. As the doppler 
measurements used a lot of spectral lines, it is not obvious to evaluate the observational 

mode mass. 

The work performed in SNS03 for the case of the Sun has been extended by I Samadi et al. 
(2005) to the case of stars lying on the main sequence for which solar-like oscillations are 
expected. They have found that the maximum in the mode amplitudes, V max , scales as 
(L/M) s where s is the slope of the scaling law, L is the luminosity and M is the mass 
of the star. The slope s was found to be significantly sensitive to the adopted eddy-time 
correlation function. As for the Sun and a Cen A , theoretical calculations assuming an 
LF fits best the maximum oscillation amplitudes detected for a set of ~ 10 stars including 
main-sequence stars as well as some red giant stars. However, these calculations rely on 
the QNA. A preliminary work tends to show that the asymmetry between the updrafts 
and the downdrafts does not change significantly between the few 3D simulations of 
main-sequence stars investigated in this work and that the CMP remains valid according 
to these 3D simulations. Hence, we expect that the difference between the effect on the 
excitation rates of the CMP model and that of the QNA remains constant for intermedi- 
ate massive stars lying on the main sequence. This however needs to be confirmed using 
more 3D simulations and extended to other domains of the HR diagram. 
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Discussion 

V.M. Canuto: How do you determine the filling factor ? 

R. Samadi: This quantity is obtained directly from the 3D simulation. 

Tim Bedding: Comment: It is very nice that you can compare your models with our 
observations. But please do not compare with amplitudes of individual peaks. There is 
a systematic effect because those peaks are selected because they happen to be highest 
during the observations. (This is true even if you average over many peaks). Instead, 
please use the amplitudes we estimated by smoothing and then noise-correcting the 
power spectrum. For a Cen A and B, this is shown in Kjeldsen et al. (2005). 

R. Samadi: Yes, I agree with you, the comparison should be done using the type of 
amplitude spectrum you mention. (Thanks to T. Bedding this is now done in the present 
proceedings). 

R.F. Stein: First, we have compared the non-adiabatic eigenfunctions from our simula- 
tions with the eigenfunctions of J. Christensen-Dalsgaard and they agree very closely, so 
I don't think that using non-adiabatic eigenfunctions is important. Second, you add the 
absolute values squared of the turbulent pressure and entropy fluctuation contributions. 
In fact, they do not always add but sometimes cancel each other. Hence, adding the 
contributions linearly would be more accurate. 

R. Samadi: We have compared the non-adiabatic eigenfunctions computed using the 
time-dependent formalism of Gabriel for convection (see Grigahcene et al. 2005) with 
adiabatic eigenfunctions and found important differences. Hence, different treatments of 
convection affect differently the eigenfunctions. Further work is thus needed on that issue. 
Concerning your remark about the turbulent pressure and entropy fluctuation, within the 
theoretical framework of our formalism, the crossing term between the Reynolds stress 
source term and the entropy source term vanishes on average. Hence, further theoretical 
work is required in order to correctly model this crossing term. 

J. Christensen-Dalsgaard: Comment: The computations of damping rates (needed 
to determine the amplitude) is probably even more uncertain than the computation of 
the energy input. For the red giant £ Hya there is a striking discrepancy between observed 
and computed damping rates. Better understanding, and better observations, are needed. 



